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2h '■ 1 Introduction 



In this paper there are studied solutions to evolutionary equations, describing wave 
patterns with compact support. Different kinds of wave patterns play key rules 
O ■ in natural processes. They occur in nonlinear transport phenomena pQ, serve as 

a channels of information transfer in animate systems [2], and very often assure 
stability of some dynamical processes [21 II]- 



^ • One of the most advanced mathematical theory dealing with wave patterns' for- 



mation and evolution is the soliton theory [5J. The origin of this theory goes back 
to Scott Russell's description of the solitary wave movement in the surface of chan- 
nel filled with water. It was the ability of the wave to move quite a long distance 
without any change of shape which stroke the imagination of the first chronicler of 
this phenomenon. In 1895 Korteveg and de Vries put forward their famous equation 

Mt + /3 M Ma; + U^xx = 0, (l) 

describing long waves' evolution on the surface of a shallow water. They also ob- 
tained the analytical solution to this equation, corresponding to the solitary wave: 

12 

u = —— sech^[a{x - Aa^t)]. (2) 

P 

Both the already mentioned report by Scott Russell as well as the model suggested 
to explain his observation did not involve a proper impact till the middle of 60- 
th of the XX century when there have been established a number of outstanding 
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features of equation ([T]) finally becoming aware as the consequences of its complete 
integrability [5]. 

In recent years ago there have been discovered another type of solitary waves 
referred to as compactons These solutions inherit main features of solitons, but 
differ from them in one point: their supports are compact. 

A big progress is actually observed in studying compactons and their properties, 
yet most papers dealing with this subject are concerned with the compactons being 
the solutions to either completely integrable equations, or those which produce a 
completely integrable ones when being reduced onto subset of a traveling wave (TW) 
solutions [IIEIE]. 

In this paper compacton-like solutions to the hydrodynamic-type model taking 
account of the effects of temporal non-locality are considered. Being of dissipative 
type, this model is obviously non-Hamiltonian. As a consequence, compactons are 
shown to exist merely for selected values of the parameters. In spite of such restric- 
tion, existence of this type of solution in significant for several reasons. Firstly, the 
mere existence of this solutions is connected with the presence of relaxing effects and 
rather cannot be manifested in any local hydrodynamic model. Secondly, solutions 
of this type manifest some attractive features and can be treated as some universal 
mechanism of the energy transfer in media with internal structure leading to the 
given type of the hydrodynamic-type modeling system. 

The structure of the paper is following. In section 2 we give a geometric insight 
into the soliton and compacton TW solutions, revealing the mechanism of appear- 
ance of generalized solutions with compact supports. In section 3 we introduce the 
modeling system and show that compacton-like solutions do exist among the set of 
TW solutions. In section 4 we perform the numerical investigations of the model- 
ing system based on the Godunov method and show that compacton-like solutions 
manifest attractive features. 

2 Solitons and compactons from the geometric 
viewpoint 

Let us discuss how the solitary wave solution to ([T]) can be obtained. Since the func- 
tion u{-) in the formula ([2]) depend on the specific combination of the independent 
variables, we can use for this purpose the ansatz u{t, x) = U{^), with ^ = x — V t. 
Inserting this ansatz into equation ^ we get, after one integration, following system 



of ODEs: 



(3) 




System (j3]) is a Hamiltonian system describing by the Hamilton function 
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Figure 1: Level curves of the Hamiltonian function of system ([3]) representing peri- 
odic solutions and limiting to them homoclinic solution 



Every solutions of can be identified with some level curve H = C. Solution 
([2]) corresponds to the value C = and is represented by the homoclinic trajec- 
tory shown in Fig. [1] (the only trajectory in the right half-plane going through the 
origin). Since the origin is an equilibrium point of system (j3]) and penetration of 
the hohoclinic loop takes an infinite "time" then the beginning of this trajectory 
corresponds to ^ = — oo while its end - to ^ = -|-oo. This assertion is equivalent to 
the statement that solution ([2]) is nonzero for all finite values of the arguments. 

Now let us discuss the geometric structure of compactons. For this purpose 
we restore to the original equation which is a nonlinear generalization to classical 
Korteveg - de Vries equation [6] : 

n, + a«™n. + /3K)_ = 0. (4) 

Like in the case of equation (P), we look for the TW solutions u{t, x) = U{C,), where 
^ = X — Vt. Inserting this ansatz into (jlj) we obtain, after one integration, the 
following dynamical system: 



dU o 

— = -nPU'W, (5) 



dT 



-VU + f/"^+^ -np W'-^W 

m + 1 



where = n P U'^j^. All the trajectories of this system are given by its first integral 



jj^+m-n _ jj5-n ^ ^ ^JJ y^,y ^ JJ ^ ^^^^^^ 



{m + 1) {5 + m — n) 5 — n 2 

Phase portrait of system ([5]) shown in Fig. [2] are similar to some extent to that 
corresponding to system (l3l). Yet the critical point U = W = of system ([5]) 
lies on the line of singular points U = 0. And this implies that modulus of the 
tangent vector field along the homoclinic trajectory is bounded from below by a 
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Figure 2: Level curves of the Hamiltonian function of system ([5]). Dashed hne 
indicates the set of singular points U = 0. 

positive constant. Consequently the homoclinic trajectory is penetrated in a finite 
time and the corresponding generalized solution to the initial system ( H]) is the 
compound of a function corresponding to the homoclinic loop (which now has a 
compact support) and zero solution corresponding to the rest point U = W = 0. 
In case when m = l, (3=1/2 and n = 2 such solution has the following analytical 
representation |6j: 



It is quite obvious that similar mechanism of creating the compacton-like solutions 
can be realized in case of non-Hamiltonian system as well, but in contrast to the 
Hamiltonian systems, the homoclinic solution is no more expected to form a one- 
parameter family like this is the case with solution In fact, in the modeling 
system described in the following section, homoclinic solution appears as a result of 
a bifurcation following the birth of the limit cycle and its further interaction with 
the nearby saddle point. 

Let us note in conclusion that we do not distinguish solutions having the compact 
supports and those which can be made so by proper change of variables. In partic- 
ular, the solutions we deal with in the following sections, are realized as compact 
perturbations evolving in a self-similar mode on the background of the stationary 
inhomogeneous solution of a system of PDEs. 

3 Relaxing hydrodynamic-type model and its qual- 
itative investigations 

We consider the following system [TO] : 





Ut+Px = 7, 
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Vt - = 0, (7) 

Here u is mass velocity, V is specific volume, p is pressure, 7 is acceleration of 
the external force, k and x/'^ ^^e squares of the equilibrium and "frozen" sound 
velocities, respectively, t is time, x is mass (Lagrangean) coordinate. The first two 
equations are convenient balance equations for momentum and mass, while the last 
one is the dynamical equation of state, taking into account relaxing properties of 
the media. 

We perform the factorization [llj of system ([7]) (or, in other words, passage to an 
ODE system describing TW solutions), using its symmetry properties summarized 
in the following statement. 

Lemma 1. System ([7]) is invariant with respect to one-parameter groups of 
transformations generated by the infinitesimal operators 



V- _d_ -y- _d_ 

^0 ~ at' 1 ~ ax' (o\ 

^ op ox oV 



Proof. Invariance with respect to one parameter groups generated by the oper- 
ators Xq, Xi is a direct consequence of the fact that system ([7j) does not depend 
explicitly on t and x. Operator X2 is the generator of scaling transformation 

u' = u, p' = e°p, V = e"" V, t' = t and x' = x. 

Invariance of system ([7]) with respect to this transformation is easily verified by 
direct substitution. 

The above symmetry generators are composed on the following combination: 



, , d d d 
[X - Xq) - — hp 7^ V ■ 



dx dp dV 



It is obvious that operator Z belongs to the Lie algebra of the symmetry group of 
(jTj). Therefore expressing the old variables in terms of four independent solutions of 
equation Z J{t, x) = 0, we gain the reduction of the initial system [II]. Solving the 
equivalent system 

dt d{xo — x) dp dV du 
T " ^{xo - x) " ^ " ^ " "0"' 
we get the following ansatz, leading to reduction: 

u = U{uj), p = Uiuj)ixo-x), V = R(uj)/(xo-x), uj = ^t + \og ^° . (9) 

Xq — X 

In fact, inserting this ansatz into the second equation of system ([7]), we get the 
quadrature 

U =^R + const. (10) 
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and the following dynamical system: 



^A{R)R' = -R [aRU - k + r^R-f] ■- 
^A{R)U' = ^ {^R {RU -n) +x{Il+l)} 



Fi, 
■■F2, 



111 



where (■)' = d (■) /dcu, A{R) = T{^Rf 



X, 



cr 



1 +re. 



In case when 7 < 0, system (ITT]) has three stationary points in the right half- 
plane. One of them, having the coordinates Rq =0, Hq = —7, lies in the vertical 
coordinate axis. Another one having the coordinates Ri = — K/7, Hi = —7 is 
the only stationary point lying in the physical parameters range. The last one one 
having the coordinates 



R2 



X 

re' 



K 



oRo 



0. 



lies on the line of singular points t{^RY — X 

As was announced earlier, we are looking for the homoclinic trajectory arising as a 
result of a limit cycle destruction. So in the first step we should assure the fulfillment 
of the Andronov-Hopf theorem statements in some stationary point. Since the only 
good candidate for this purpose is the point A {Ri, IIi), we put the origin into this 
point by making the following change of the coordinates X = R — Ri, Y = 11 — Hi 
which gives us the system 



X 
Y 



— K 

?2 



-Rja 



X 
Y 



Hi 
H2 



;i2) 



where 



Hi = - [UiX^ + 2aRiXY + aX^Y) , 
H2 = {UiX^ + 2RiXY + X^Y) . 
A necessary condition for the limit cycle appearance read as follows [T2] : 



spM = ^ {^Rif + x^ 



detM >0 ^ = <A(i?i) > 0. 



(13) 
(14) 



where M is the linearization matrix of system f[T^ . The inequality f[T^ will be 
fulfilled if ^ < and the coordinate Ri lies inside the set (0, ^/x/ (t^^))- Note that 
another option, i.e. when ^ > and A > is forbidden from physical reason 
On view of that, the critical value of ^ is expressed by the formula 



2Rl 



(15) 



Remark. Note that as a by-product of inequalities inequalities ([T3|) . f[T4j) we get 
the following relations: 

- 1 < < 0. (16) 
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To accomplish the study of the Andronov-Hopf bifurcation, we are going to cal- 
culate the real part of the first Floquet index Ci [12j. For this purpose we use the 
transformation 

enabling to pass from system (fT2l) to the canonical one having the following anti- 
diagonal linearization matrix 

Mij = Q{62i6ij - 6ii62j). 

For this representation formulae from [121 ^^^e directly applied and, using them 
we obtain the expression: 

IGRln^Re Ci = -K {3K^ + {^Rif (3-^r) -k R^f (6 + rO} ■ 

Using ( fT3l) . we get after some algebraic manipulation the following formula: 

Re C, = Y^Y^ {2^^^^ (^^i)' - (^'^0' - 3 (XO'} • 

Since for ^ = < expression in braces is negative, the following statement is 
true: 

Lemma 2. \i Ri < R2 then in vicinity of the critical value ^ = E,cr given by the 
formula ( JTSl) a stable limit cycle appears in system ( fTTl) . 

We've formulated conditions assuring the appearance of periodic orbit in proxim- 
ity of stationary point v4(i?i. Hi) yet in order that the required homoclinic bifurca- 
tion would ever take place, another condition should be fulfilled, namely that on the 
same restrictions upon the parameters critical point B{R2, II2) is a saddle. Besides, 
it is necessary to pose the conditions on the parameters assuring that the stationary 
point B{R2, II2) lies in the first quadrant of the phase plane for otherwise corre- 
sponding stationary solution which is needed to to compose the compacton would 
not have the physical interpretation. Below we formulate the statement addressing 
both of these questions. 

Lemma 3. Stationary point B{R2, II2) is a saddle lying in the first quadrant 
for any ^ > C,cr if the following inequalities hold: 

-T^,rR2 < Rl < R2. (18) 



Proof. First we are going to show that the eigenvalues A 1^2 of the system's f[TT|) 
Jacobi matrix 



M 



d{F,,F2) 
d {R, n) 



eHa-2)+2'y Rija^l)] _xo_ 
rr ' T 



(19) 
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are real and have different signs. Since the eigenvalues of M are expressed by the 
formula 



sp M ± y sp M - 4 det M 
Ai,2 = ^ 

it is sufficient to show that 



detM < 0. (20) 



In fact, we have 



detM 



X 



[k (fr-2) + 27i?2 (a-1)] 



r T 




) 



2x^7 (^1-^2) < 



To finish the proof, we must show that stationary point B{R2, II2) lies in the 
first quadrant. This is equivalent to the statement that 



Carrying the first term into the RHS and dividing the inequality obtained by 7 < 0, 
we get the inequality —T^crRi < Ri- The latter implies inequalities — r^i?2 < 
Ri < R2 which are valid for any C, > ^cr- And this ends the proof. 

Numerical studies of system's flTT]) behavior reveal the following changes of regimes 
(cf. Fig. [3]). When ^ < ^cr, A{Ri,^i) is a stable focus; above the critical value a 
stable limit cycle softly appears. Its radius grows with the growth of parameter 
^ until it gains the second critical value ^cr2 > ^cr for which the homoclinic loop 
appears in place of the periodic trajectory. The homoclinic trajectory is based upon 
the stationary point B{R2,Il2) lying on the hne of singular points A(i?) = so it 
corresponds to the generalized compacton-like solution to system ([7]). We obtain 
this solution sewing up the TW solution corresponding to homoclinic loop with 
stationary inhomogeneous solution 



corresponding to critical point B{R2,Il2)- So, strictly speaking it is different from 
the "true" compacton, which is defined as a solution with compact support. Note, 
that we can pass to the compactly supported function by the following change of 
variables: 

Trit, x) = pit, x) — 112 (a^o ~ a;), u{t, x) = V{t, x) — R2/ {xq — x). 

4 Numerical investigations of system (|7|) 

4.1 Construction and verification of the numerical scheme. 

We construct the numerical scheme basing on the S.K.Godunov method |I5j. Since 
the inhomogeneous terms appearing in ([7]) destroy the scaling invariance, we look. 



K, — T ^crl R2 > 0- 



u = 0, p = Il2{xo-x), V = R2/{xo - x) 



(21) 
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Figure 3: Changes of phase portrait of system (fTTj) : (a) A{Ri, Hi) is the stable focus; 
(b) A{Ri, Hi) is surrounded by the stable limit cycle; (c) A{Ri, Hi) is surrounded 
by the homoclinic loop; (d) A{Ri, Hi) is the unstable focus; 



in accordance with common practice [TB], for the solution of the Riemann problem 
{ViUiPi) at a; < (V2, at x > to corresponding homogeneous system 



Ut + Px = 0, 

Vt-u, = 0, (22) 
X_ 



Using the acoustic approximation, we find the functions U, P in the sector —Ct < 

x<C%C= ^xli-rViY 

U="-l±^^Eigi, (23) 

p ^ P1+P2 ^ Ui - U2 

2 2 

where Vq = ^^^^^ . Expression for the function V is omitted since it does not take 
part in the construction of the scheme on this step. 

At some additional assumption the Riemann problem can be solved without re- 
sorting to the acoustic approximation. Let us assume that 

As it easily seen, this relation is the particular integral of the third equation of 
system f l22]) . Employing this formula, we can write down the first two equations as 
the following closed system: 

1.4) (M.o, ....A-JX-^'r)^ <-) 
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Solving the eigenvalue problem det \\A — = 0, we find that the characteristic 
velocities satisfy the equation 

Now we look for the Riemann invariants in the form of infinite series 

oo 

It is not diffiecult to verify by direct inspection that the following relations hold: 

D^V = ± Cloo^^ V = u,± ClooV, = Q±, (26) 

D^u = ±ClooQ±- (27) 

Using (12^ . we find the recurrency 

and finally obtain the expression for Riemann invariants: 

r^=AoVexp {^u/^^). (28) 

So under the assumption that p = system (1221) can be rewritten in the following 
form: 

D±r^ = 0. (29) 



Using and fl2^ we get the solution of Riemann problem in the sector ^/x/(T^^t < 
x<Vx/FW)t: 



U = In Z, (30) 

P = P2 + V^C2[Zexp i-U2/^/^) - 1], 

where = v^/^^ = Cloo(V^«), ^ = 1,2, 



Z = (E+ v^)/(2C2Vx7^), 

E = exp (wa/ \/x/t~){Pi ~ P2 + \/x/r{C2 - Ci)}, 

Q = E'^ AxlrCiC2 exp [(m + u^)/^/^]. 

Note that ( 130|1 is reduced to ( l23l) when |pi — P2I << I7 |wi — U2\ « 1. 
The difference scheme for ([7j) takes the following form: 



[u7 - ur')Ax - (K+1/2 - K-i/2)Ai = -lAtAx, 



(V^r - Vr')Ax + («IVi/2 - <_V2) A^ = 0' 
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where (m"_i/2' Pl-1/2) Pl+1/2) are solutions of Riemann problems (V^li, m"-!) Vi-i] 

(V;", <, p^) (V;", M^, p^), {V^l-^, M^Vi' ^r+i). correspondingly, 

i 

k is equal to either n or n + 1. The choice k = n leads to the explicit Godunov 
scheme 

vr^' = vr + - (31) 

The scheme (1311) was tested on invariant TV solutions of the following form: 

u = U{uj), p = P{uj), V = V{tu), uj = x-Dt. (32) 

Inserting fl32|) into first two equations of system ([7]), one can obtain the following 
first integrals: 

U = u, + D{Vi-V), 



P = p^ + D\V,-V), 

where Vi = lim^^^ao ^(^)- Let us assume in addition that Ui = 0, while pi = n/Vi. 
With this assumption constants {ui, pi, Vi) satisfy the initial system. 
Inserting U and P into the third equation of system ([7]) we get: 

where Ctoo = a/x/t, S = pi + D^Vi. Equation (j34l) has three critical points: 

V = Vo = 0, V = Vi, V = V2 = k/{ViD^). 

If the inequality V2 = k/{ViD'^) < Vi holds and the line x/r - (DV)'^ = is 
outside the interval then (V2, Vi) then constants 

(DViY - 

U-00 =U2 = — > 0, p_oo =P2 = K./V2, I/-00 = V2, 

deliver the second stationary solution to the initial system and solution to fl34|) 
corresponds to a smooth compressive wave connecting these two stationary solutions. 

Results of numerical solving the Cauchy problem based on the Godunov scheme 
( 13T|) are shown in Fig. |H As the Cauchy data we took the smooth self-similar 
solution obtained by numerical solving equation and employment of the first 
integrals (!33ll . So we see that the numerical scheme quite well describes the self- 
similar evolution of the initial data. 
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Figure 4: Temporal evolution of Caucliy data defined by solutions of equation 
and the first integrals (l33l) . Following vaues of the parameters were chosen during 
the numerical simulation: k = 0.5, x = 0.25, r = 0.1, Vi = 0.5, D = 3.1. 

4.2 Numerical investigations of the temporal evolution and 
attractive features of compactons. 

Below we present the results of numerical solving of the Cauchy problem for system 
([7]). In numerical experiments we used the values of the parameters taken in accor- 
dance with the preliminary results of qualitative investigations and corresponding 
to the homoclinic loop appearance in system (ITT]) . As the Cauchy data we got the 
generalized solution describing the compacton and obtained by the preliminary solv- 
ing of system ffTTl) and employment of formulae ffTOj) . fl2T]) . Results of the numerical 
simulation are shown in Fig. [51 It is seen that compacton evolves for a long time in 
a stable self-similar mode. 

Additionally the numerical experiments revealed that the wave packs being cre- 
ated by sufficiently wide family of initial data tend under certain conditions to the 
compacton solution. Following family of the initial perturbations have been consid- 
ered in the numerical experiments: 



Here a, /, pi, w, h are parameters of the perturbation defined on the background of 
the inhomogeneous stationary solution fl2Tl) . Note that parameter / defines the width 
of the initial perturbation. Varying broadly parameters of the initial perturbation, 
we observed in numerical experiments that, when fixing e.g. value of Z, it is possible 
to fit in many ways the rest of parameters such that one of the wave packs created by 
the perturbation (namely that one which runs "downwards" towards the direction 
of diminishing pressure) in the long run approaches compacton solution. Whether 




Po{xo — x) when x G (0, a) U (a + /, Xq) 
'i){xo — x) + w{x — a) + h when x G (a 
u = 0, V = n/p. 




(35) 
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Figure 5: Numerical solution of system ([7]) in case when the invariant homoclinic 
solution is taken as the Cauchy data 



the wave pack would approach the compacton solution or not occurs to depend on 
that part of energy of the initial perturbation which is carried out "downwards". 
Assuming that the energy is divided between two wave packs created more or less 
in half, we can use for the rough estimation of convergency the total energy of the 
initial perturbation, consisting of the internal energy Eint and the potential energy 



E. 



pot- 



E = Eint + E, 



-/pot 



pot. 



dx, 



where Eint, and Spot are local densities of the corresponding terms 

Function Spot is connected with forces acting in system by means of the evident 
relation 

_ ldp_ _ depot 

p dXe dXe 

where Xg is physical (Eulerian) coordinate connected with the mass Lagrangean 
coordinate x as follows: 

X = dxp. 



Prom this we extract the expression 



E.. 



V 



dp 

dXe 



7 dx'e 



Cl 



V~^dx, 



where Q is the support of initial perturbation. Employing in the above integral the 
relation V dp/dxe = dp/dx, we obtain the following formula: 



pot 



a + P 



1 + k 
k 



ln(l + A;) - 1 



where = [P(a +/)-P(a)]/P(a), P{z) = az + (3, a=w-{po+Pi), (3 
{pi +Pq)xq - aw. 
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Figure 6: Perturbations of stationary invariant solutions of system ([71) (left) and TW 
soluations created by these perturbations (right) on the background of the invariant 
compacton-like solution (dashed) 
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Figure 7: Evolution of the wave patterns created by the local perturbations which 
do not satisfy the energy criterium 



For X = 1-5, K = 10, 7 = —0.04, r = 0.07, Xq = 120 convergency was observed 
when Efot was close to 45 (see Figures below) . 

Function Emt is obtained from the second low of thermodynamics written for the 
adiabatic case: {deint/dV)s — —p — —k/V. From this we get 



'int 



To obtain the energy of perturbation itself, we should subtract from this value the 
energy density of stationary inhomogeneous solution c — «;ln Vq, so finally we get 

Eint = J i^int - el^^)V~^dxe = n j \n.Vo/Vdxi. 



Using the forula (3. 7. 20), we finally obtain 



Eint = S / In — -rr H lu 

Po(a + 1) « 



1 + 



al 



a, 



P2 



1 



P2I 

Po{a) 



where Po{z) = po{xo - z). 

Numerical experimenting shows that the energy norm serves as sufficiently good 
criterium of convergency. At x = 1.5, k, — 10, ^ — —0.04, r = 0.07, Xq — 120 
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convergency was observed when E G (43, 47). Patterns of evolution of the wave per- 
turbations is shown in Figs [61 For comparison we also show the temporal evolution 
of the wave packs created by the perturbations for which E ^ (43, 47) (Figs. [7]). 

Thus there is observed some correlation between the energy of initial perturbation 
and convergency of the wave packs created to the compacton solution. 

5 Conclusions and discussion. 

In this work we have discussed the origin of generalized TW solutions called com- 
pactons and have shown the existence of such solutions within the hydrodynamic- 
type model of relaxing media. The main results concerning this subject can be 
summarized as follows: 

• A family of TW solutions to ([7]) includes a compacton in case when an 
external force is present (more precisely, when 7 < 0.) 

• Compacton solution to system (JTj) occurs merely at selected values of the 
parameters: for fixed k, 7 and x there is a unique compacton-like so- 
lution, corresponding to the value ^ = ^cr2- 

Qualitative an numerical analysis of the corresponding ODE system describing 
the TW solutions to initial system served us as a starting point in numerical inves- 
tigations of compactons, based on the Godunov method. Numerical investigations 
reveal that compacton encountering in this particular model form a stable wave 
pattern evolving in a self-similar mode. It is also obtained a numerical evidence 
of attracting features of this structure: a wide class of initial perturbations creates 
wave packs tending to compacton. Convergency only weakly depend on the shape 
of initial perturbation and is mainly caused by fulfillment of the energy criterion. 
This criterium is far from being perfect. In fact, it is not sensible on the form of ini- 
tial perturbation, which, in turn, influences the part of the the total energy getting 
away by the wave pack moving " downwards" . Besides, employment of the Godunov 
scheme does not enable to obtain more strict quantitative measure of convergency. 
But in spite of these discrepancies the effect of convergency is evidently observed 
and this will be the topic of our further study to develop more strict criteria of 
convergency as well as trying to realize whether the compacton solution serves as 
true or intermediate [TTl [IB] asymptotics. 
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